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Abstract. We compute the third order term in a generalization of the Strong 
Szego Limit Theorem for a zeroth order pseud odiffe rential operator (PsDO) on 
a Zoll manifold of an arbitrary dimension. In |G02 , the second order term was 



computed by V. Guillemin and K. Okikiolu. In the present paper, an important 
role is played by a cert ain c ombinatorial identity which we call the generalized 
Hunt-Dyson formula [ Gli3 1 . This identity is a different form of the renowned 



Bohnenblust-Spitzer combinatorial theorem which is related to the maximum 
of a random walk with i.i.d. steps on the real line. A corollary of our main 
result is a fourth order Szego type asymptotics for a zeroth order PsDO on the 
unit circle, which in matrix terms gives a fourth order asymptotic formula for 
the determinant of the truncated sum Pn{T\ + T2D)Pn of a Toeplitz matrix 
Ti with the product of another Toeplitz matrix T2 and a diagonal matrix D 
of the form diag(0, li ■ ' ■ )• Here Pn = diag{l, ■ ■ • , 1, 0, • ■ ■ ), n ones. 



1. Introduction 

The main motivation for this work was to find an exphcit formula for a "Szego- 
regularized" determinant of a zeroth order pseudodifferential operator (PsDO) on 
a Zoll manifold introduced in |G01, after (3)] and [02|, see Remark 1.7. Our 
main result, Theorem 1.3, is valid for any dimension 0? G N. In the case d = 2, 



Theorem 1.3 gives such a formula. 



1.1. Notations and main results. Let M = §^ be the unit circle R/27rZ. Denote 
by P„, n e N, the orthogonal projection from L^(§^) to the subspace spanned 
by {e'^^^jifeKn. For a function / € L^(S^) denote its fcth Fourier coefficient by 



fk ■— Jq^ fix)e If, fc e Z. Let b{x) be a positive function on such that 

Sfeez 1^1 l(loS^)fcP < Denote by B the operator of multipHcation by h acting in 
L^(§^). The matrix representation of the operator B in the basis {e^^^^k^i is the 



Toeplitz matrix {bj- 
states that 



c)j,feez- The classical Strong Szego Limit Theorem (SSLT) 



00 



TrlogP„BP„ = TrP„(logB)P„ + > A: (log 6)^ (log5)_^ + o(l), n ^ cx). 



fe=i 



Here TrlogP„PP„ = logdetP„PP„ and TrP„(logP)P„ = (2n+l) /p'"log6(x) 
It h as b een shown by H. Widom that the remainder is 0{n^°°) if h{x) £ C°°(S^), 
see [Wll. 
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The main result of this paper is Theorem 1.2, in which we find a third order 
generahzation of the SSLT for a zeroth order pseudodifferential operator (PsDO) 
_B on a ZoU manifold M of an arbitrary dimension d S N. 

Recall that M is called a Zoll manifold |G01| if it is compact, closed and such 
that the geodesic flow on M is simply periodic with period 27r. The unit circle 
and the standard sphere of any dimension are Zoll manifolds. A second order 
generalization of the SSLT for a Zoll manifold M of any dimension has been obtained 
by V. Guillemin and K. Okikiolu |G01, G02], see also an important preceding 
work @ by K. Okikiolu for M = S^md S^The proofs in |0l|, |G01| , |gU2| use a 
combinatorial identity due to G. A. Hunt and F. J. Dyson and proceed in the spirit 
of the combinatorial proof of the classical SSLT by M. Kac See also [G03, p2| 
where the combinatorial approach and the Hunt-Dyson formula (HD) are used in a 
different setting to obtain a second order generalization of the SSLT for a manifold 
with the set of closed geodes ies o f measure zero in the u nit cot angent bundle. 

In the proof of Theorem L3 we use the method of | G02 |. A central role in 
our proof is played by a certain combinatorial identity which generalizes the Hunt- 
Dyson formula mentioned above to an arbitrary natural powe r. W e call this identity 
the generalized Hunt-Dyson formula (gHD), see Theorem 8^ and Gi3 |. After 



having discovered and proved the gHD we realized that it is related to another 
combinatorial theorem, which has a long history. This theorem is a result due to 
H. F. Bohnenblust that appeared in an article by F. Spitzer on random walks |^ 
Theorem 2.2], and is now commonly known as the Bohnenblust- Spitzer theorem 
(BSt). A major application of, and motivation for the BSt, is the computation of 
the characteristic function of the maximum of a random walk with independent 
identically distributed (i.i.d.) steps carried out in Note that the expectation of 
such a maximum was computed earlier in |Q with the help of the usual HD. 

Let M be a Zoll manifold of dimension d G N. Let ^'™(M), m e Z, denote the 
space of classical PsDO's of order m on M. Recall that for a given G G ^'™(M), 
its principal symbol am{G) and subprincipal symbol sub(G') are well-defined on 
T*M. Let A denote the Laplace-Beltrami operator on M. It is known |DG| that 
there exists a constant a £ R such that the spectrum of A lie s in bands around 
the points fc + j, fc e N. Moreover, it has been shown in [ |CdV | that there exists 
j4_i G ^'"^(Af) such that [A, = and the spectrum of the operator 



(1.1) 



A := A_ 



is N. Let Pn, n e N, denote the projection from L'^{M) onto the subspace spanned 
by the eigenfunctions of A corresponding to the eigenvalues 1, 2, • • • ,n. Let dxd^ 
be the standard measure on S*M :— {(x,^) : ai{A){x,£^) — 1} divided by {2TTy. 
Following | G02 we will assume that (Ti{A){x,£,) — ai{A){x,—£_) for all (a;,^) £ 
T*M. In I Gil, Chapter 1] this is not assumed which leads to more complicated 
expressions. Let 6*(a;, ^) denote the shift of the point {x, £,) £ S*M by t units along 
geodesic flow. For any function / e C°°{S*M) introduce the fcth Fourier coefficient 
along the closed geodesic of length 27r starting at a given point (x, ^) 



(1.2) 



dt 



e~^^'f{Q\x,i))—, keZ 
Zn 



The simplest form of our result is for the case of M = §^ with f{z) = logz. 
Note that f{z) is analytic in a disk of radius 1 about the point z = 1. In the proofs 
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in Section ^, we require that the function f{z) is analytic on a disk the radius of 
which depends on a certain norm of the operator B E 5'°(M). For our purposes 
the following norm is convenient 

\\\B\\U := ||i?|| + \\[A, [A,B]]\\ + ||V"^'^(™)ao(B)||oo 
^^■^^ +||V3sub(B)|U+ / |vV_d(A2-'^B)|da;de, 

where d ~ dimM, || • || is the operator norm in L'^{M), V includes both x- and 
^-derivatives in local coordinates on T*M, and ||g||oo ■— inax(^x,i)eS'M\g{x,£,)\- 
(The integral in is well-defined being a Guillemin-Wodzicki residue, see ( ^.81 ) 
below.) 

Theorem 1.1. Let M ~ and Pn he the projection on the linear span o/{e*'^'^}|fc|<„ 
Let B G ^^{M) and assume that <Jq{B) is strictly positive, and that the symbolic 
norm |||/ — _B|||i is sufficiently small. Then \ogB G ^^{M) and the following holds 
as n oo, 
(1.4) 



TrlogP„BP„ = TrP„(logB)F„ + - / V fc ((To(logS))^ {ao{\ogB))_^ dxd^ 
f Y.k (aoTbTi?)),. (sub(k^S))„, dxd^ + o(-. 



1 1 

n ' 2 



In (1.4) the argument (x,^) G 5**^/ is omitted for brevity and for each [x,^) G 
S*M the Fourier coefficient is understood in the sense of (1.2). 

We need to fix more notation to formulate the result for M — §^ and an arbitrary 
analytic f{z). Let Ai denote the set of all analytic functions on C with no constant 
term 



Ai {/(z) : f{z) = ^ ^ '^}- 



In I LR£ ] the authors introduce a linear map W2 from Ai to the space of continuous 
functions from to C, defined by 

(1.5) W^mx^x^) := \ r r tM^lMd^^d^,. 

^ Jo Jo ?i ~ ?2 

Let j G N. We will call a linear map V from Ai to the space of continuous functions 
— !■ C a linear j-map. The linearity means that for arbitrary /, (/ G Ai, a, /5 G C, 
xi, ■ ■ ■ ,Xj e C 

V[af + (3g]{xi, - ■ ■ ,Xj) ^ aV[f]{xi, ■ ■ ■ ,Xj) + (3V[g]{xi, ■ ■ ■ ,Xj). 



A 2-map U, which is equivalent to W2, was earlier constructed by H. Widom |W6 
W7|. We will need also a 2-map W2 whose action on an arbitrary / G -4i is 
prescribed by 

(1-6) W2[f\{xi,X2) := - dii. 

2 Jo ?i - ^2 

For an arbitrary B G '^^{M), let us write 60 := cro{B) and 6sub := sub(_B), 
&o(x,^) := 6o(e*(a;,$)), hl^^{x,S,) ^sub(6*(a;, ^)), and omit the argument (x,^) G 
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S*M. It is convenient to introduce the notations 

(.271 /.2ir 



•J S* AI ^ -j^ "^0 "^0 



and 



Theorem 1.2. Let M = and P„ he the projection on the linear span o/{e*'^^}|fc|<„. 
Let B e *°(M) and f e Ai. Then f{B) G ^°{M) and the following holds as 
n CO, 



1 



T3,sub[/](S)+0 



1 



Observe that Theorem |l.l| follows from Theorem |l.2| setting f(t) = logz and 
noting following H. Widom and [LRS, p2| that 

(1.7) W2[log](xi,a;2) = logxi loga;2, 

and also that for < ti,t2 < 2tt 

WMibl\bl')bll^ - i^ogbl'MlJbl') = iMlog B)Y- (sub(log 

In the higher dimensional case two additional contributions to the third Szego 
term now arise. First, define a 3-map W3 such that for any / € Ai 



W3[f]{xi,X2,X3) := X3 



(1.8) 



Jo 

/fe) 



ei(a-2;3)(6-6) 

/(^3) 



6(6 - 2^3) (6 - 6) 2:3(6 - 2:3) (6 - X3) 

and introduce the notation 



p2Tr /.27r 



T3,o[/](B):=(d-l) / dxd^ 

f2(k' + il + a/2)k) r I 
^ -'o Jo 

- E rrr^"^'^"^'^'^*^"*^^^3[/](&^6^6^^)fif^f^ 

fc,i=i Jo Jo n n 77 



see Remark |1.2| . Secondly, we introduce notations needed to describe the contribu- 
tion of the Poisson brackets {&q, < t, s < 2tt, to the third Szego term. Define 
for each j £ N a linear j-map $j such that for any m £ N, m > j , 

(1.9) $,[z"](:ri,...,x,):= ^ 



and ^^[z™] := for TO=l,---,j — 1. This together with the linearity defines 
uniquely on Ai, see Section \i\ a nd (7.2) for an expression for <&_,■ acting on an 
arbitrary / g Ai . We can write ( |l.9| ) in terms of the complete symmetric function of 
degree m — j evaluated at the point (6, • • • ■ ■ ■), see (7.1). Now for x,y £R 
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denote —{x)- := min(0, x) and M2{x, y) := min(0, x, x+y). For arbitrary j, /c, / G N 
and any xi, X2, Mi> • • • , fJ-j, i^ir ■ ■ ,J^k, Pi,-- - , Pi & set 

Ctf\xi,X2,IJ,l,--- ,P'j) 

:= M2{xi,X2) - (xi + X2-M2(xi,X2)- 

(1.10) - (>.i)_ - ( (>.i)+ - (Ml)- (W)- 

- ((Mi)+ + --- + (Aij)+ + >*2)_)_ 

- (Mi)- (Mj)- - ((m)+ H 1- (Mj)+ + ^2(>^l,>f2))_, 

and 

^j\fc(^i'^2,Mir • • t^ir- - ,f/c) 

:= -(^fi)- - ((^fi)+ - (/xi)- inj)- 

- ((/il)+ + ••• + (/ij)+ +X2 - {Ul)- (^'fc)-)_)_ 

- (/il)_ ifXj)- - ( (/il)+ + • • • + + XI - (X2)- 

(1.11) -((X2)+- (:.!)_ (^.)-)_)^ 

- (mi)- (Mj)- - ( (mi)+ H H (/Uj)+ + XI 

- (!/l)_ (z/fc)_ - h (t'fe)+ +X2)_)_, 

and 

^jfM(^i'^2,Mi,--- ,fJ'j^^ir-- ,i^k,Pi,--- ,Pl) 

■■= -(Ml)- (Mj)- - ( (Mi)+ + • • • + (Mj)+ 

(1.12) + XI - (z/i)_ (!/fc)_ - ( + ■■■ + {Uk) + 

+ X2 - (pi)- (pi)_ )_ 

For any f G Ai denote by Tj[f]{z) its Taylor polynomial of degree j G N about the 
point t = 0. For any f G Ai and B G ^^{M) introduce the notation 



T 



3,Poi[/](5) := /^^^ (a(i)[/](&o) + A(2)[/](5o) + A(3)[/](6o)) dxd^ 



where 



11 

A(^^[/](M ^ Et? E 0«(xi,X2,Mi,--- ,Mi) 

j = l '<i+x'2+(i*iH \-Hj=0 

(1.13) xiT'"""""""'f«'«"'^^ 

,.277 ^2ir 

oi(AtiriH hMj-rj) 



Jo Jo 



X [.-^(/(.) - T2[/](.))] (6-, • . . , |i . . . ^, 
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and 



(1.14) 





+ ^l + -- + >'.-=0 



Jo 

27r <.27r 



dui du2 



dr\ drj dsi 



gi(AiiriH hMj-fj+i^isiH hi^kSk) ^ . . . 3 

27r 27r 27r 



27r 



and 



' "0 ^ 



2i ^ 7! fc! l\ 

j,k,l=l 



E 

+X2 +^^1 H l-Atj + 

+ ,^j + ... + ,^j,+pjH hPi = 

/■2-K f2ix 



(1.15) 



X / / e 

'0 Jo 

X 



i(^lUl+^2ll2) r /,"2-l ^^^2 



27r 



dr\ drj ds\ ds^ dti dti 



Jo 

1(^1 i-iH hAiji"j+i^isiH hi^fcSfc+pitiH Vpiti) 



X e 



' "0 ' "0 ' 



Now the result for any dimension and an arbitrary / G Ai . Denote 
T3[/](i3) T3,o[/](S) + T3,.ub[/](B) + T3,Poi[/](i?). 

Theorem 1.3. Let M be a Zoll manifold of dimension d e N. Let A be defined 
by (|L1]). Assume that (7i{A){x,^) = ai{A){x,~-^) for all (x,^) G T*M. Let B e 
^"{M) and f £ Ai. Then f{B) G ^"(Af) and the following holds as n — > 00, 



(1.16) 



TTf{PnBPn)^TrPnfiB)Pn 



T2[/](B) 

Tsifm + Oin"-^). 



Most of the paper is devoted to the proof of Theorem |1.3| . Theorem |1.2| follows 
from Theorem p^.3| in view of the following. For d = 1, T3.0 vanishes, and also 
Tr3,poi vanishes, because all the Poisson brackets vanish in this case (for each of the 
two cotangent directions the angle does not change and ao{B) is homogeneous of 
degree in ^). Note that for the function f{z) — logz, the terms involving W2 and 
W2 in Theorem 1.3 for any d G N take a simpler form as in Theorem 1.1. We make 



several remarks, and then state a corollary to Theorem 1.1 and 1.3 which gives an 
explicit formula for logdetP„i?P„, as n ^ cxd. For d = 1, 2, that formula gives an 
expression for a possible generalized determinant of the operator B G 4'°(Af), see 
Remark 1.7 below. 
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Remark 1.1. The existence of a full expansion of the type (1.16) for f{z 
TO G N, f(z) ~ logz, has been proven in 



two coefficients were given in |G01, p02 
f{z) = log 2;, and for general f e Ai, in LRS|] 



See Remark 2.2 in subsection 2.2 
for a discussion of the formulas for further coefficients. 



G01 |. Explicit expressions for the first 
in the case when f(z) — z'", m G N, 



Remark 1.2. The formula (LS) for M^3[/] has a structure similar to the coefficient 
in the third asymptotic term in aSzego type expansion for convolution operators 
obtained by R. Rocc aforte in |RR|. A certain combinatorial identity known as 



Spitzer's formula (see DyMcK , Section 3.7]), which is a ver sion of the Boh ncnblust- 
Spitz cr co mbinatorial theorem, is also used in the proof in | RR |. (See also Bax , Ro| 



and Gi3 , Remark 1.6]). A second order Szeg o type expansion for convolution 
operators was established by H. W idom in |W9| with the help of the usual Hunt- 
Dyson combinatorial formula (8.4). A full asymptotic expansion for convolution 



operators is obtained in |W4] 



Remark 1.3. If / G Ai is a polynomial then each of the sums over j,k,l in 
T3,poi[/](-B) reduces to a finite sum for any ZoU manifold M and any given B G 
^°(M), because $m[z"] = identically for m > n. The term T3,Poi[/](-B) van- 
ishes for any / G .4i in the case of c? = 1, and also in the following case: Assume 
/(^) = E^i C2fc-iz2'=-\ let M S'^, d G N, d > 2, and let B G ^'°(M) have a 
principal symbol bo G C°°{S'^) which is independent of ^, and in addition is an odd 
function on Thus bo{x) — —bo{xa) for all pairs of antipodal points 
Then the Poisson brackets in T3_Poi[/](i3) are odd with respect to ^, and vanish 
after the integration over S*M. (We refer to [W3| for an auxiliary calculation of 
the needed Poisson bracket on S^, which can be easily modified for any W^, d > 3.) 
Under these assumptions the coefficient of the first asymptotic term in ( 1.16| ) of 
order n'^, which is a part of Tr Pnf{B)Pn, also vanishes, but not the coefficients of 
the second and the third asymptotic terms of orders n'^~^ and n"^"^, respectively. 



Remark 1.4. In order to define the third Szego asymptotic term in (1.16) the con- 
dition / G C'*(/) is necessary, where / is the closed set of values of a-o{B). If the 
term T3_Poi[/](i3) is abse nt (se e Remark 1^) then the condition / G C^{I) is nec- 
essary. It was shown in [ LRS that if B is self-adjoint and /" G L°°{I) then the 
second order Szego formula holds, see also |LS|. It would be interesting to know if 
the condition of essential boundedness of the fourth derivative of / is sufficient for 
Szego asymptotic formula with three terms to hold. Another question concerns the 
best possible norm in place of (pT^). 

Remark 1.5. We have discovered and proved the generalized Hunt -Dyson combi- 
natorial formula (Theorem ^.2[ ) being unaware of the Bohnenblust-Spitzer theorem 
(Theorem 8.1). A derivation of the gHD from the BSt and vice versa can be found 
in [ pi3[ ], see |Gil| , Chapter 2] for an independent proof of the gHD. Also the impor- 
tance of the BSt in the theory of the maximum of a random walk with real i.i.d. 
steps and related results are discussed in | Gi3 , Remark 1.5 and 1.6]. 



1.2. Explicit asymptotic formulas for logdetP„i3P, 
logz. Theorem 



1.1 



and 1 . 



as n - 

give an expression for logdetP„BP„ 



00. Let f{z) = 
= Tr log P„BP„ 

as a sum of TrP„(logP)P„ and two lower order corrections, as n — s- 00. Propo- 
sition 1.4 below gives, for an arbitrary G G ^'''(M), an auxiliary asymptotic ex- 
pansion for Tr PnGPn, as n — > 00. For dimension d = 1,2, we need the constant 
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coefficient in this expansion, which is more compHcated than the other ones, see 
the proof of Proposition lA in Section || for de tail s. Let Ri{G), I — 0, 1, 2, • • • , 
be the Guilleniin-Wodzicki residues as given in (2.8) below. For d — 1,2 make an 
additional assumption 



(1.17) 



1=0 



C{G) :=^(Tr(^fcG)-5]fc'^-i-'i?,(G) 



under which the following sums are absolutely convergent 

(1.18) " ' 

fc=l 

see Section || for the proof. 

Let 7 denote the Euler constant and C, the Riemann zeta function. 

Proposition 1.4. Let M be a Zoll manifold of dimension c? G N. Let Pn he as 

above and assume that G € 4'°(Af). For d — 1,2, assume in addition that (1.17) 
holds, and let C{G) he defined by ( |1.18 ). Then the following holds as n ^ oo, 
(i) for d = l, 

TrPnGPn ^n-Ro{G)+ logn • i?i(G) 



+ (g(G) +7i?i(G) 

(a) for d — 2, 

TrPnGPn^n^ ^ 



E 

1=2 



ai)Ri{G) 



1 

n 



i?i(G)-i?2(G) +0 



2^o(G) 



RiiG) 

O 

(G(G)+7i?2(G) 



log71-i?2(G) 



Y,ai)Ri+i{G) 

1=2 



o 



0' 



(Hi) for d > 3, 

Tr PnGPn = 

'd-1 

■ n" 



12 



1 
' d 

Ro{G) 



Ro{G) 



^'-'■{iRoiG) 



d- 1 



i?i(G) 



■i?i(G) 



d-2 



i?2(G) +logn-i?<j(G) +0(n 



d-3\ 



Remark 1.6. The coefficients of n'^, n*^" -^ and log n for d> 2 can be found in | G02 
after Lemm a 0. 2]. F rom Proposition 1.4 with G = logi? we see that Tr P„ (log i3)P„ 
in Theorem 1.1 and 1.3 contributes to the leading asymptotic term of order n'', and 
also to all lower order terms of order , j — d — 1, • • • , 1, 0, — 1, • • • , and to the 
logarithmic term logn, as rt — + oo. In the classical SSLT the situation is much 
simpler: logi? is just the Toeplitz matrix of the operator of multiplication by log 6, 

and so TrP„(logB)P„ = (2n + l)(bg6)o. 

Now we are ready to state two corollaries. 

Corollary 1.5. Let B £ ^''^(S^) have a strictly positive principal symbol and with 
a sufficiently small symbolic norm |||/ — P|||i. Assume also that (1.17) holds. Then 
the following holds as n ^ oo, 



log det PnBPn = ci ■ n + ci^g ■ log n + cq 
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where the coefficients are the sums of the corresponding coefficients from Theo 
rem 



1.1 and Proposition 1.4(i) for G = \ogB. 
Assume further that aa{B) and sub(-B) do not depend on the direction of S^, that 
is ao{B){x,£,) = bQ{x) and anh{B){x,^) = 6sub(a;) \^\^^ , for (a;,^) S S*Ei^ . Assume 
also that 6_2 ~ 0. Then the following holds as n ^ oo, 
(1.19) 



logdetP„BP„ = n • 2 
+ log n ■ 2 



I ^ogbo{x) — 



27r 



bn{x 

OO 

^fc(log &o)fc(log&o) 



C(logB)+7i?i(logB) 



. fc=i 



E 

1=2 



C{l)Ri{\ogB) 



^A:(log6o)fc(6sub/&o) 



, k=l 



p2TT 




f bsnh{x) 


)1 


dx \ 


lo 


bo{x) 


\ bo{x) 




2^ j 



01 4 



where C{logB) is given by (1.18) 



The proof of ( pj| ) is an exercise m the calculus of PsDO's together with Propo- 



sition L4(i), and is left to the reader. In so me simple cases, for instance for 
bsuh{x) = ±i&o(a^)j the left-hand side in ( 1.19 ) can be computed explicitly. The 
coefficients of n, log n, and ^ on the right in ( 1.19 ) in these cases are as expected, 
see also Remark |l.9| below. 

Corollary 1.6. Let M be a Zoll manifold of dimension d > 2. Assume that P„ 
and A are as in Theorem l.S. Let B € ^''^(Af) have a strictly positive principal 

\L — B\\\i. For d = 2, assume in 
~* oo, 



symbol and a sufficiently small symbolic norm 
addition (1.17). Then the following holds, as n 



(1.20) log det P„SP„ -C^' 



c[^l-\ogn + 0(n 



d-3\ 



wher e th e coefficients are the sums of the corresponding coefficients from Theo- 
rem LS for f(z) — log z and Proposition \l.^ (ii) or (Hi) for G — logP. // one 
counts the logarithmic term, this expansion is fourth order for d — 2,3 and third 
order for d > A. 



Remark 1.7. The coefHcients G^f^ and C^"\, deN, have been found in pOl|,|G02| 



id) 



The most interesting coefficient in (1.20) is the constant one, since one can think 
of expCg'''' as a regularized determinant of P e 5'''(A/), see [ pOl , after (3)] and 
(0|. The sum 

OO 

7Pd(log B) + Y, Cmi+d-i{\og B) 



1=2 



will for all d G N be a part of C^f^ . For d = 1, Corollary gives a full expression 
for C^^\ For d = 2, Corollary [O gives a full expression for Cq^\ which is quite 
lengthy. 
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Remark 1.8. Let us compare the result of Corollary 1.5 with a generalization of 
SSLT to the case of B being an operator of multiplication by a function h{x) having 
discontinuities which is due to H. Widom and E. Basor. In this case log b{x) also has 

discontinuities, and so the series X^fcez 1^1 l(log6)fcP diverges logarithmically. The 
following third order asymptotic formula holds for the operator of multiplication 
by a piecewise function h{x) 

(1.21) log det P„i?P„ = ai • n + 02 • log n + as + o(l), n ^ oo, 



where ai as in (|1.19| ), the coefficient 02 has been computed by H. Widom in [|W6| , 



and the constant term 03 has been found by E. Basor in M. Note that the matrix 



B in (1.21) is still Toeplitz, the logarithmic order of the subleading term being due 
to a slower decay of the Fourier coefficients of b{x) . In our case the matrix of the 
operator B G \l/°(§^) is not Toeplitz (see Remark |l.9| below), and the logn term 
comes from the contribution of sub(i3). 



It would be interesting to find a compact formula for the constant term in (1.19) 
We mention that the constant 03 in ( |l.21[ ) found in |^ has a form similar to the one 



in (1.19). It contains a "finite" term and and an infinite series of certain integrals 
multiplied by the values of the Riemann zeta function at the points 3, 5, • • • . Inter- 
estingly, an "invariant" form of that series has been found in | |W^ . It is written as 
a single integral involving the function 

*W ^ ^ogV{x). 



This gives the hope that a similar formula can be found for the constant in ( 1.19 ) 



Remark 1.9. The matrix interpretation of Corollary L5 is as follo ws. Assume 
for simplicity that B e 5'°(S^) is as in the second part of Corollary L5, that is 
CTo(S)(a;,0 = bo{x) and sub(B)(x,^) = hnh{x) |Cr\ for aU (x,C) e S'*§^ Assume 
also that 6_2 ~ ^-3 = • • • = 0. Let Bq and i?sub be the operators of multiplication 
by 60 and 5sub, respectively. Let D be the linear operator in L^(§^) such that 

|0, fc = 0. 

Note that this is not a differential, but rather a smoothing operator of order 
— 1. There is known a correspondence between the classical PsDO's on the cir- 
cle and their discrete counterparts, see @ for details. By that correspondence, 
the zeroth order PsDO B we started with equals i?o + B^uhD. Introduce two 
Toeplitz matrices, Bq {(&o)j_fc}j,feez and Bsub {(&sub)j_fc}j,fcez- Set also 
D := diag(- • • , i, i, 1, 0, 1, i, i, • • • ). Then the matrix representation of Bo+BsuhD 

is Bo + Bsuh ■ D. Finally, set P„ = diag(- • • , 0, 1, • • • , 1, 0, • • • ) (2n 1 ones). Then 
Corollary |l.5| gives a fourth order asymptotics for the determinant of the truncated 
matrix P„ • {Bq + B^ub ■ D) ■ P„. 

Now we can reformulate the question of finding the constant term in ( 1.19] ) in 



purely matrix terms. Drop the hats and the dots for brevity. Let Ci be a Toeplitz 
matrix that corresponds to the operator of multiplication by 6sub/&o, and let the 
matrix D be as above. Clearly, the matrices Ci and D do not commute. Assume 
that the matrix log(/ — CiD) is well-defined. The question is to compute the 
constant coefficient in TrP„log(J — CiD)Pn, or which is the same, the constant 
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coefRcient in 



Tr Pnlog{I - D^/'^ClD'^''')Pr 



As we have noticed in Remark 1.6, this question is trivial for a Toephtz matrix T 
in place of D^/^CiD^/^ . 

1.3. A related result for the maximum of a random walk. Let us explain 
how we can use the j-maps j £ N, to write the bivariate characteristic function 
of the maximum of a random walk and its position at a smaller time. There is a 
lack of symmetricity in this problem, and this case is not considered in |^]. Let 
Xi, X2, • • • be independent real valued random variables which assume real values 
and have the same distribution density (l){X). Assume for simplicity that is 
Schwartz class. (The result below holds for much more general (p, e.g., for discrete 
random variables, if understood in the sense of distributions.) In our case the 
characteristic function 

/•OO 

E{e'*^} := / e''^(j){X)dX 



is well-defined and invertible as a Fourier transform. Let Sp := Xi + • • • + Xp be 
the position at time p G N of the random walk starting ofl^ at the origin. For any 
g S N introduce the random variable 



p+q 



(0, S*!, • ■ ■ , S'j), S") 



p+l> • 



which is the length of the maximal excursion to the right of during the time 
interval 0, 1, • • • ,p + q. Note that the time p + q is strictly larger than p. Let us 
introduce a non- negative valued function similar to the functions ( 1.10| ), (1.11) and 
( 1.12| ). For arbitrary j, fc G N and yi, - ■■ , yj, zi, - ■ ■ , g R we set 



(1.22) 

Denote 

0(77) 



%fc(yi,- • 

+ ( 



• • • ,Zfc) := (2/1)+ H h {yj)+ 

(yi)_ {y,)^ + {zi)+ + --- + izk)+)^. 



e-^^yc^{y)dy, .^-^^/(r;)] (2^) 



Theorem 1.7. Let (f>, Xp, Sp, and Tp^g, p,q £ N, be as above. Let $j, j G N, be 



as in ( |l.9| ). Then for \a\,\b\ < 1 the following holds 

00 



p,q=l 



(1.23) 



00 ^ 

j\ki 



X $1 



J.a{yx^ |-l/j)+i/5f2j,j.(ai,--- ,yj-Z\,--- ,Zk) 



bt ■ 



as 



1-bt 



1 - 

(0(Ci), 



(0(771),- •• ,0(?7j)) 



,0(a)) 



dyi-- ■ dyj dzi ■ ■ ■ dzk 



In this formula, the j-maps $j and $fc act on the function of s and t, respectively. 
The arguments of these actions are the values of (j) at the corresponding points. 

The proof for fixed p, g G N is carried out analogously to the computations in 
Section |[ see [Glil, Section 1.10] for details. 
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1.4. Organization of t he art icle. The paper is organized as follows. In Section || 
we recall the method of G02 |, estimate the remainder after the third Szego term, 
and justify the passage from the set of polynomials to an analytic function /. 
After that we deal with an arbitrary monomial f{z) — z™, m g N. An operator 
5'°(Af), >?i, • • • , >c„i £ Z, arises in that computation (see (2.3) below 



for the definition of the Fourier coefficient i?^, >f € Z). In Section^, we compute the 
contribution of (Jq^B^^-^ ■ ■ ■ i?>c„) , to > 2, to higher order Szego terms. The resulting 
expression involves only cto(-B). In Section ^, we calculate the contribution of the 
symmetric part of suh(^B^^ ■ ■ ■ B^^^, m > 2, to the third Szego term. It depends 
on (7o(-B) and sub(i?). In Section ^ we compute the contribution of the non- 
symmetric part of snh^B^-^ ■ ■ ■ B^^ ) , m > 2, to the third Szego term, which involves 
Poisson brackets of the principal symbol of B shifted along by geodesic flow. This 
contribution depends only on CTo(i?). In Section |^, a proof of Proposition 1.4 is 
given. In Section |^, we find an expression for j G N, and also for Waf/], for 

an arbitrary / e Ai . In Section ^, the aux i liary co mbin atorial background is given 
T he r esult s of this paper (Theorem 



1.1 



and |8.2|, Proposition 1.4 and Corol 



lary L5 and 1^) were announced in [ Gi2 |, where we also gave an outline of the 
proofs. 

Acknowledgments. This paper is a modified part of the author's Ph.D. the- 



sis I Gil I . I would like to express a deep gratitude to my thesis adviser Ari Laptev, 
for suggesting the problem and his constant attention to the work. 

I would like to thank Percy Deift for giving a reference |RS| from which I learned 
about 1^ and the Bohnenblust-Spitzer combinatorial theorem, shortly after having 
discovered and proved the generalized Hunt-Dyson formula. 

I would also like to thank Kurt Johansson for suggesting a derivation of the 
formula for f E Ai, in Section 0, via generating functions. 

Finally, I would like to thank the Department of Mathematics of the Royal 
Institute of Technology (KTH) , Stockholm, for a generous financial support during 
the whole period of my graduate studies. 

2. The full third Szego term and the remainder estimate 



2.1. The abstract scheme from [G02]. We start by e xpanding the analytic 
function f{z) in a power series about and proving ( 1.16 ) for f{z) — z™ for an 
arbitrary to G N. After that we justify the passag e from the set of polynomials 
to analytic functions. Let us recall the method of |G02|. Let TTfc, A: S N, be the 
projection on the fcth eigenspace of the operator A and set tt^ := for < 0. Then 
P„ = J2k=i ^fc foi" n£N, and we set P„ := 0, n < 0. For an arbitrary B £ ^'''(Af) 
and t eR introduce the operator 

By Egorov's theorem, i?* G ^*'(Af), and also 

(2.1) ao{B*){x,0^'Jo{B){eH^,0), 

where 0* stands for the shift by t units along the geodesic flow. Note that because 
sub (A) = const the following also holds Q 

(2.2) snh{B'){x,0 = sub(B)(e*(x,e)). 

Because spec(A) = N, the operator i?* is periodic in t with period 27r. Therefore 
we can introduce the Fourier expansion 

B = Exez where B^ e^°(M),}€& Z, 
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is defined by 

(2.3) = E S TTfe = — / e^-* e-^'^Be^'^ dt. 

This together with (|2l|) and ^) imphes 
Lemma 2.1. i^or any S e *°(Af) anrf e Z 



and 



sub(B.)(:E,0 - e'"* sub(i?)(e*(a;,e)) ^. 

For TO €E N and , • • • , >c„-i G Z introduce the notation 

(2.4) Mm(Jc) min(0, m:i + >i2, ■ ■ ■ , H 1- >*m)- 

Now using the remarkable commutation relation 

(2.5) B^Pn - Pn+^B^, n e N, ^ G Z, 
we move all the projectors to the left in the expression 

(P BP — \^ PR P B P ■ ■ ■ P R P 



obtaining PnB^"'Pn plus another term. This implies for all n G N 
Tr(P„BP„)'"-TrP„B"/'„ 

= — Tr ((Pn ~ PnPn+>ci ' ' ' Pn+niA \->Cm)P>^l ' ' ' P >'m) 

(2.6) >^i+---+>c„.=o 



= - E E Tr(7r„+fcS^i •••B>^„). 

>tlH hXn,=0 j = A/„(^) + l 



Remark 2.1. The relation ( p.5[ ) follows readily after writing both sides using the 
definition ( p.3[ ). We note that the Fourier coefhcients B^t are PsDO's even if B 
is a multiplication operator. Also even in the simplest case M — S^, B ~ b{x), 
B^ 7^ 6;j^e*^^. This explains the non-symmetricity of ( p.5| ) with respect to positive 
and negative 

Next, for any G G ^"^(M), M being a ZoU manifold, there exists a full asymp- 
totic expansion for Tr(7rfcG'), as fc ^ oo, see Lemma 2.2 below. This result is 
due to Y. Colin de Verdiere |CdV |. The coefficients in that expansion are certain 
Guillemin-Wodzicki residues. Recall that for any compact closed manifold M of 
dimension d G N the Guillemin-Wodzicki residue of a pseudodifferential operator 
G G \I'™(A'f) of order to G Z is defined by 

(2.7) Res{G) :=[ cr^4G){x,0 dxd^ 

(recall that we have included {2tt)^'^ in the notation dxd^). For an arbitrary G G 
^^i^iM) denote 

(2.8) Ri{G) ■.= Rcs{A-'^+'G), l = 0,l,2,---. 
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Lemma 2.2. Let M he a Zoll manifold of dimension d £ N. Assume G G \I'°(M). 
Then for any N = Q,1^2, ■ ■ ■ , there exists Cn < oo such that 



(2.9) 



N 



Tr(^fcG)-^fc'^-i-'i?KG) 



1=0 



fc e N. 



See |CdV| | and ]G02| , Appendix] for the proof of (|j|). We need the expressions 
for the first two residues. 



Lemma 2.3. Let A be defined by (1.1). For any G e ^'"(M) 



(2.10) 

and 

(2.11) 



Ro{G) 



CTo(G) dxd^ 



S'M 



Ri{G) = / 

J S'M 



{d - l)(a/4) (7o(G) + sub(G) ) dxd^. 



Proof. The equahty ( 2.10| ) follows easily from the computation rules for PsDO's 
and the fact that (Ji{A) = 1 on S*M. Next, the equality 



(2.12) 



i?i(G) = / svih{A-'^+^G)dxd£^. 

J S'M 



follows from the definition of the subprincipal symbol (see, e.g., [DG|), and the fact 
that the integral over S*M of a Poisson bracket equals zero. (We could also refer 
to Proposition 29.1.2 in [H64|, differentiate with respect to the spectral parameter, 
make the same rema rk con cerning the integral involving the Poisson bracket.) Now 
( ^.ll| ) follo ws fr om ( ^.12| ) by the computation rules concerning the subprincipal 
symbol: by pCl (1.4)] 



sub(A 



-d+l, 



'G) = sub(A-''+^) ao(G) + a^d+iiA-''+') sub(G) 
+ (2z)-i{a_,+i(A-'^+i),ao(G)}, 
where {•, •} denotes the Poisson bracket. Now by | |DG| , (1.3)] 

sub(A-'*+i) = (-d+l)(ai(A))"'^ sub(A), 
and we use the definition ( |1 . 1| ) , and that sub(— A) = to complete the proof. □ 



2.2. Beginning of the proof of Theorem |l.3| . We do not assume that the 
Fourier expansion B = "^^i^^B^ has only a finite number of terms. Denote 



(mi,-- - ,Xm), 



Bjt :— B^-^ ■ ■ ■ 



Below c{B) and C{B) will denote various constants which depend only on B (and 
not on to). We will indicate which seminorni of B enters a certain C{B) when 
necessary. For any M > denote 

(2.13) 

\^\<M 

Introduce for m G N, M > the set of indices 

g'"(M) :={(><i,... ,,^,„) gZ" : \Hi\<M,l^l,--- ,m}. 



LOWER ORDER TERMS IN SZEGO TYPE LIMIT THEOREMS ON ZOLL MANIFOLDS 15 



We split the right-hand side of (2J) into the following three sums. The first sum 



IS 



(2.14) 
the second 
(2.15) 

and the third 
(2.16) 



E 



-eQ'"("/(2™)) , = M^(>c) + l 



E E Tr(7r„+j Bjr) 

(n + jf-^ Res{A-'^Bj;) - (n + jf-"^ R%s{A^-'^B^) 



E 



E Tr(7r„+j 



The contributions to the second and third Szego term come from (2.14), the ex- 
pressions ( ^.15 ) and ( 2.16| ) having a lower order, as n ^ cxi. In the rest of this 
section we estimate the remainder after the third Szego term. It is very important 
that the symbol of the operator B of any order is smooth. We will repeatedly refer 
to subsection 2.7 where the most technical part of the computation is carried out. 



2.3. Compu tatio n of the sum ( 2.14 ). Let us single out the terms involving n'^ ^ 
and n''^^ in ( 2.14 ). Then the latter splits into the sum of 



(2.17) 



E E 

«eQ"(>«/(2™)) i=M„(>?) + l 



( (d - 1) • j • R(is{A-'^B^) + R&%{A^-'^B^ 



and the remainder 



r d-3 



(2.18) 



>C-l-\ [--Wt71=0 



E E Y.^H\)3'-^-''-R^<A-^'B^) 



d-3 



d-2 



Taking into account Lemma 2.3 for G = B— and the formulas 



J2 1 = -Mm( 





j=M^il<) + l 



(M„(3<))2-f M,„(3<) 
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we conclude that ( ^.17[ ) equals 



E 



-eQ'"(n/(2m)) 



S'M 



(2.19) 



d-2 



d- 1 



(Af„(>?)) +(1 + -)M„(>?) 



S'M 



+ Mm{'>i) I snh{B^)dxdC 

Js* M 

In Section ^, the expression 

(2.20) V (M„ (>?))" / (7o(G^i---G^„)dxd^ 

for n G N and any G G 4'°(M) is computed (we set G := _B("/(2™))). The impo rtant 
point here is that both the domain of summation and the second factor i n (2.20 ) 
are symmetric, and so we can symmetrize the first factor. The integral in (2.20) is 
symmetric with respect to >fi , • • • , >Cjn because already the integrand is symmetric 

o'o(G^i • • • G^^) = o'o(G^i) • • • cro(G^,„). 



The gHD for n = 1, 2 (Theorem B^) is needed in this computation. 
In Section and pi the expression 



V M™(>?) / suh{G^,---G^Jdxd^ 



>flH |-J<m=0 

for any G G ^'"(M) is computed (we set G := b("/(2™))). Note that 



sub(G^, • • • G^J = sub(G^J n MB>.J 



k=l 



(2.21) 



p=i 



2i 



Y {cTo(G.J,co(G.J} n ao(G.J 



l<fc<i<m 



p=i 



where the first sum is symmetric with respect to >fi , • • • , >Cm , and the second one is 
generally speaking not (not even after the integration over S*AI). Because of that 
circumstance the original method of symmctrization |Q ^, |0l| , G02| , G03, p2| 
fails, and we have to modify it. It turns out that each of the m(m — l)/2 terms 
in the second sum in ( ^.21 ) possesses a partial symmetry. For instance, if 1 < 
r + 2 < s and s + 2 < m then in the expression 

o'o(Gj<i • • • G^^) {(To(G^^^i), cro(G>f^^J} o'o(G^^^2 ' ' ' C^s) o'o(G>c^^2 • • • G^^) 

the indices >fi , ■ ■ • , >^r can be permuted without changing the resulting expression, 
and the same holds for the groups of indices >Cr+2, ■ ■ ■ , Xs and >fs+2, • • • , (and 
we are even allowed to interchange the indices within the three groups). However 
we can neither interchange any index from any of the three groups with >Cr+i, Xs+i^ 
nor interchange the latter two. Here the original form of the BSt (Theorem ^.l[ ) is 
needed. 
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Remark 2.2. In view of ( p.6[ ) and Lemma 2.2, the computation of all terms in 
the Szego asymptotics is reduced to the evaluation of the following expression for 
different n e N and / = 0, 1, 2, • • • , 



(2.22) 



d+l 



■ ■ ■ B^, 



The gHD (Theorem 8.2) makes the computation possible for any n e N, provided 
that the second factor in (2.22) is symmetric with respect to xi,--- ,Xm- The 
problem is that this is the case generally speaking only for I ~ 0. The case of Z = 1 
and n = 1 is dealt with in Section ^ and ||. 

We would like to mention that in the work | RR | by R. Roccaforte an interesting 
way to rewrite expressions of the type ( |2.22D with non-symmetric (matrix-valued) 
second factor has been suggested. The idea is to consider all the cases when the 
minimum is attained on the sum 



xi 



J = 1, • ■ 



1, 



and make certain changes of summation indices. More precisely, for an operator 
B e ^''^(Af) written as a Fourier series B — X^xgz following |RR| we set 



B, 



B'+^ := {B * Bl) jeN 
B'+^ -.^{B^^B) jeN, 



where * denotes the discrete convolution on the Fourier series side. The order of 
the operators being convolved is important. Then ( 2.22] ) can be rewritten as 



(2.23) 



^ {-{x)-T Res 



A 



-d+l 



B 



Note that the summation in (2.23) is over a single variable k. This reminds the 
expression which appears in the usual second Szego term 1^, |G02|] Also the 
formulas from ||W2}] have a similar structure. However the formula (|2.23|) for /(z) — 
z™ is not very explicit, and we cou ld not write a reasonable formula for an arbitrary 
analytic f{z) starting from ( 2.23| ) (nor was it done in | RR | ) . 

Bringing together the results of Section 3, 4, and 5, we conclude that ( p. 19 ) 
equals 



(2.24) n'^-i • T2[2"](s("/(2m))) _^ „<i-2 . -[^\^z"'^^^B 

This in turn equals 

nd-i . T2[z'"](B) + n"-'' ■ T3[z'»](B) 
plus an error whose absolute value does not exceed 



(n/(2m))^ 



n 



(2.25) 



T2[z'"](5^"^^^'"") - T2[z"](B) 



,d-2 



T3[z'"](B("/(2™))) - T3[z"](B) 



n e N. 



We need the following auxiliary statement which is proved in subsection 



2.7. 
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Lemma 2.4. For I = 2,3, an arbitrary B E ^'''(M) and any L E N there exist two 
constants ci.l{B) and Ci{B) such that for m E N, m > 2 



Tz[z™](i3(^^))-T;[z'"](i3)^ 

< ci^l{B) ■ ■ {CiiB)r ■ M-\ 



M > 1. 



It follows from the proof of Lemma that Ci{B), I = 2,3, involve s || V ^'7o(i?)Hoo 
and II sub(i?)||oo, but does not depend on L. Applying Lemma 2A with AI = 
n/{2m) and I — 2, L = 2, respectively, I — 3, L — 1 to the first, respectively, second 
term in (2.25), we estimate the latter by 



c-m''-c(B)-(C(B))", neN. 



where C{B), 1 = 2,3, depends on \\\7^<Jo{B)\\^ and || sub(B)||oo. 

We finish this subsection by showing that the absolute value of the remainder 
( 2.1S| ) can be estimated, for n G N, by 



n'^-^ ■ c{B) ■ m ■ {C{B)y 



For instance, the absolute value of each of the d ~ 2 terms in the first sum in the 
square bracket in ( 2.18[ ) does not exceed 

J2 iMn^m" I Res{A-"B,, 



n^-^-Cd 



(2.26) 



E 



{1 + l^ilf ■ ■ ■ (l + \^„-,\f 



(^1,-, 



X / dxd^\ao(B^^)\ 

J S'M 

< n'^-^ -ciB) • m - (C(B))" 



\MB>cJ\ 
n eN. 



The last inequality is due to the fast decay of the maximum over S*M of |g o(-B> c)|, 
as |>!:| ^ oo (smoothness of cro{ B) an d Egorov's theorem), see subsection 2/7 for 
details. We remark also that in ( 2.26| ) C{B) depends on || V'^+^o-o(-B)||oo- In the 
same way one shows that the absol ute value of each of the d—2 terms in the seco nd 
sum in the square bracket in ( 2.1§| ) is estimated by the right-hand side of ( 2.26 ). 

2.4. Computation of the sum ( ^.161 ). For any G Z™ \ Q"'{n/{2m)) the abso- 
lute value of at least one of the components of >f is > n/{2m). Redenote it by >ri. 
We also note that j > — n in (2.16) (otherwise iTn+j = 0). Then the absolute value 
of the sum (2.16) does not exceed 

E 







m J2 

j— — n ^GS 



'":|«ll>n/(2m) 
l + -- + «m=0 



(2.27) 



Tr(7r„+jB^j - --B^^) 
E Cd{n + jf-^\\B^ 



\B, 



«i + -- + «„,=o 



< n ■ m ■ Cd 



E 



Eii^^ 



7n — 1 



|>ci|>n/(2m) 
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where the factor Cd{n + j)'*^^ is an estimate of the niultiphcity of the eigenvalue 
n + j and || • || is the operator norm in L?'{M). By ( 2.39 ) below, there is a constant 
c{B) such that 

-3 



E 

\xi\>n/(2m) 



< <B) — 
' 2m 



n 



n e N. 



Set C{B) := E^ezM^W- % {M3 below, C{B) < \\B\\ + (7rV3)||[A \A,B]]\\. It 



follows now from ( p.27j ) that the absolute value of the sum ( 2.16 ) is estimated, for 

71 e N, by 



CdciB)-m^-{C{B)y 



2.5. Computation of the sum ( 2.15| ). For an arbitrary G £ ^'''(Af) and any 
/I € N denote 



(2.28) 



,-d+3 



Tr(^^ G) - fi"-' Res(^-'^G) - fi"-' Res(A^-''G) 



,d-2 



We wish to prove that there exist t wo co nstants c{B) and G{B) so that the limit, 
as n ^ oo, of the absolute value of ( ^.15|) d ivided by n'^~^ is < c{B) ■ m ■ {C{B))™. 
Let us estimate the absolute value of ( 2.15| ) divided by n'^^^ by 



E 







(2.29) 



>^lH \-J€jn=0 



n + j 



< 



E 



rf-3 





E 



Tn+jB-^ 



Tn+jB^ , 



^1 H |-^TTT=0 



where we have used that n/2 < n + j < n for >5 e (5™(n/(2m)). We need to 
interchange the lim„_^oo to the right-hand side in (2.29) and the sums over infinitely 
growing sets. To justify that we find a summable over and j majorant and then 
refer to the Lebesgue dominated convergence. We use the principle of uniform 
boundedness to present a summable majorant. 

Introduce a Banach space X = 5'°(Af)/5'~'^(Af) with a norm || • \\x given by 



X 



\\G\\+ / {\a_d{A-'G)\ + \a_d{A'-''G)\ 



+ \cr_d{A^-'^G)\)dxd^, GeX, 



where t he in tegrands are well-defined by the definition of the Guillemin-Wodzicki 
residue ( ^.7[ ). Then it is clear from the definition ( 2.28|) that each : X ^ C is 
a linear and bounded functional, because the multiplicity of each eigenvalue /l( is 
finite and 



Res(^'-'^G)| < / I a-diA^^^G) \ dxd^, 1^0,1, 2. 



Also by Lemma | 
(2.30) 



for all G e X there exists the limit 
lim T^G = Res(A2-'*G). 



Therefore for all G G X one has sup^^p^ I^^Cj < oo. By the principle of uniform 
boundedness, there exists a constant N < oo such that 

sup|T^G| < Ar||G||jf, GeX. 
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Wc finish the construction of the majorant by showing that it is summable. The 
right-hand side in ( 2.29| ) is estimated by 





B^W^ < N 



E 



|M™(3?))| \\By, 



«e«'"(»i/(2m)) 



I I ) II \\^ < OO, 



<N {\>,,\ + - 

where the last inequaHty is due to the smoothness of the symbol of B of any order, 
see subsection 2.7. 

The passage to the limit under the infinite sums in the right-hand side of ( p. 29 ) 
is therefore justified. We have 





lim 

n — >C30 



(2.31) 



E E 

-eQ™(i/(2m)) j=Af„(>f)+l 

XI H ^^m=0 



Tn+jB-^ I 



E 



I -i [->i-T7i=0 



Elim 
n— *oo 

i=A/,„(l<) + l 



7ri+j ( i? 



?("/(2m)) 



where the notation ( 2.13| ) has been employed. Due to t he in finite smoothness of 
the full symbol cf{B) and Egorov's theorem, and also by (|2.30| ) 



lim 



Tn 



Q(n/(2m)) 



hm \ TnBj;\ 



Res {A^^'^B, 



dx d£, 



Therefore the right-hand side of ( 2.31 ) does not exceed 

V {\^i\ + --- + \><m\) I a^4A^-''B^,---B^ 

< c{B)-m-{C{B)r, 

again due to the smoothness of the symbol of any order of B, see subsection | 

We notice that C{B) depends on /^.^^^ \V^c7_d.{A^^'^B)\ dxd^, || V^cto(-B)||oo and 
II V3 sub(B)||oo. 



2.6. End of the proof of Theorem 1.3. We have shown that there exist two 
constants c{B) and C{B) (the latter constant is small together with the norm (1.3)) 
such that for all m € N, to > 2, one has 



(2.32) 



Tr((P„SP„)'") - Tr(P„i?™P„) - n"-' ■ T2[z"](P) 



,d-2 



•Tain (5) 



<n'^-^-c{B)-m''-{C{B))"', n ^ oo. 



Take now an arbitrary function f{z) — X]m=i ^m-z'" € which is analytic on a 
neighborhood of {z : |2;| < C{B)}. Because the trace is a linear operation and the 
functionals T/, I — 2, 3, arc linear in the first argument, we can write 

Tr( /(P„PP„) ) - Tr(P„/(P)P„) - n'-' ■ T2[/(z)](P) - n^-^ ■ T,[f{z)]{B) 

OO 

J2 Cm [Tr( {PnBPr^r ) - Tr(P„P™P„) 



m— 1 
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Now we estimate the absolute value of the latter sum by the sum of the absolute 
values of its terms and employ ( 2.32 ). We obtain then 

Tr( fiPnBPn) ) - Tl{Pnf{B)Pn) - u"'^ ■ T 2[f]{B) - u"'^ ■ T3 [/] (B) 
00 

the series on the right is convergent due to the analyticity of / on a neighborhood 
of {z : \z\ < C{B)}. 



2.7. Proof of Lemma |2.4| . We have to prove that for an arbitrary B e '^'^{M), 
1 — 2,3, and any X e N there exist two constants ci^l{B) and Ci{B) so that for 
TO e N, m > 2 

< ci,l{B) ■ ■ {Ci{B)r ■ M-"-, M ^ oo. 
We consider the case I — 2 first. By the definition of T2 

T2[z'"](B(*^')-T2[z"](S)| 

V M^(k) I cjo{B^'^ ~ Bjr) dxd^ 



(2.33) 



We have 



(2.34) 



B.M'P-B..)Bi'^^---Bil) 



Let us estimate the integral of the first summand, the other m — 1 are estimated 
analogously. Note that Af„j(x) < (l + |><i|)---(l + \ xm\). Also because the sum 
of the indices in (2.33) is zero we have 

- Xl Xm\ <\X1\^ 1- \Xm\ 

and hence 

Af„,(>?)| < (l + |>^2|)'---(l + |><,n|)'. 

(We do not want >c\ to appear on the right-hand side in ( 2.3E1 ) below.) Introduce 
the notation for G £ ^'°(Af) 

||cro(G)||oo ^ max |cto(G)(x,0|. 

Hence 

M™(>?) / MiBi^P - B^jBi^'P ■ ■ ■ B\^) dxdi 

>ciH |->cm=0 JS'M 

(2.35) <\S*M\ (l + l><2l)'---(l + k™l)' 



)\\J\MBi'P)\U---\\a,{BZ')\ 



For any G £ \I'"(A/) we write for brevity 50 ■— o'o(G), gsub ■— sub(G), introduce 
the notation 



(2.36) 



(AdA)(G) := [A,G] e ^%M) 
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and state the following important fact (see | G02 , Lemma 1.3]). 
Lemma 2.5. For any G G \E'°(Af) and any S Z 
(2.37) vG,^[A,G],, 

where the right-hand side is the vth Fourier coefficient of [A, G] . Also for any L G N 
there exist finite cl{G) and Gl{G) such that 



and 
(2.38) 

and 

and 



\>t\>M 



||sub(G,)|| <CLiy' 



■ CX), 

M oo. 



J2 lisub(G^)||oo <ClM-^, 

\>t\>M 

and for the operator norm in i^(M) 

(2.39) \\G^\\ <GlM-^, M^oo. 

\x\>M 



Proof. The equation (2.37) follows readily from the definition (2.36) and the fact 
that TTp is the projector on the corresponding to the eigenvalue /x € N eigenspace 
of the operator A. Due to the infinite smoothness of ao{G) we can repeat the 
operation (2.37) any finite number of times, and consequently the norm ||(t(G'i/)||oo 
decays rapidly, as oo. By the fast decay of ||cro(Gj^)|joo, as |>f| oo, for any 

L g N there is a constant Cl(G) such that for aU {x, ^) € S*M and <t <2tt 



= \ Y 'To(G.)(e*(x,o 

< Y lko(G^)||oo <Gl(G)-M-^, M^«3 



\>c\>M 



which together with Lemma proves (|2.38| ). The result for the su bprin cipal 
symbo l also follows from its infinite smoothness and Lemma 2A. For ( 2.39 ), see 
G02| , Lemma 1.3]. □ 

For any G G *°(M) we now define inductively (AdA)^(G), N e N, and note 
that 

(Ad A)^(G,) = ((Ad A)^(G))„ J/ e Z, TV e N. 
For = 1 this follows readily from ( ^.36 ), for A^ > 2 we proceed by induction. 

Lemma 2.6. For any H e *°(A/) and £ Z 

lko(i?.)IU < lko(i?)||oo 

and 

\\suh{H^)\\oo < ||sub(iJ)||oo 



Proof. Follows from Egorov's theorem {2A) and (|2.2| ) 



□ 
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Lemma 2.7. For any G G 'I'°(M) 

(2.40) J2 lko(G-)IL ^ lko(Go)|loo + y ||fTo((AdA)2(G))| 



^ II sub(G,)||^ < II sub(Go)||oo + y II sub((Ad^)2(G))| 
and for the operator norm in L'^{M) 

(2.41) E 11^-11 ^ l|Go|| + y ||(Ad^)2(G)||. 



Proof. Applying (2.37) twice we get 

which imphes ( |2.40| ) because X^jy^^o '^"^ = '^^Z^- For ( |2.4l| ), see |G02| , Lemma 1.3]. 

□ 

Now we estimate ( ^.35 ) by 

Elko((i?(^^)-s).JL 



( 5: (1 + |>.2|)^|ko(i?if ))||oo) • • • ( E (1 + M?h,{Bil))\\. 



By ( 2.38| ) for any L e N there is C2.l{B) such that, as M ^ oo, 

E lko((s(*'^-i?).jL- E lko(s.j|L<c2,L(i?)M-^ 

>tiGZ \>ci\>M 

Also 

5] (1 + k2|)'|ko(B£f ^)||oo = E (1 + 2^21 + \^2\')\WoiBi'P)\\. 

and for example 

\>C2\<M 

= ^ ||ao((Ad^)2(i?))^J|U 



(2.42) 



|3<2|<M 

< ^ |ko((AdA)2(i3))^J||^=:G(i3), 
3<2ez 



where G(i3) is estimated by || V'*cro(i?)||oo, in view of Lemma 2/7. We have proved 
that 

T2[zn(s(*')) - T2[z™](S) I < C2,l(B) • m • (^2(5))™ • M"^, M > 1. 
where Ci{B^ depends on || V*(To(-B)||oo- 
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Now the proo f for T3. In view of ( 2.19 ), the functional T3, which is the coefficient 
of n^^'^ in ( 2.24), c ontains three parts. The first one is up to a constant as the right- 
hand side of ( ^.33| ), and is estimated in the same way as above. The second part 
is 



V (M„(3?))2 / ao{B!i'^ - B^) dxd^ 

-n JS*M 



>tlH h^m=0 

and the only difference with the above argument is that here we have (1 + |>£';|)^, 
I — 2, - ■ ■ , m in the estimates, and consequently 



CiB) :=5]||ao(((AdA)4(B)), 



76, 



can be estimated in terms of |jV (To(^)IIoo, by Lemma 2.7. The absolute value of 
the third part of T3[z™] (S^^)) - T3[z'"](B) is estimated by 



V Af™(>?) / sub {b!^''^ - B^) dxd£. 



>^iH 1->^™=0 



We use (2.34) and estimate for instance the integral of the first of the m summands 



sub ( [B 



(M) 



B.,)Bl'P ■ ■ ■ BZ' ) = snhiB^^^> B^Jao{Bli-> ■■■B 



AM) ^ _ 



AM) 



AM)] 



(2.43) + MBi'P - S.J sub ) . . . Bii^)) 

+ lr{ao{Bi^^-B.J,MBi^y--Bi^)} 
Denote G := (Ad A)'^{B). By Lemma for any i e N we obtain 

5^ M„(55) / snh{Bi:p^B^,)ao{Bi'P---Bi'^))dxdi 

<CL-m-{C{G))"' ■ M-^, M>1. 

We note that C{G{B)) can be estimated in terms of || sub(i?)||oo and || V'*cro(i?)||oo' 
Now for the Poisson bracket term in (2.43). 

Lemma 2.8. For any G G *"(M) and L e N there exist cl{G) and Cl{G) such 
that 



and 
and 

and 



00 



Ci(G) :=^||9.ao(G^)||oo <oo 
G2(G) :=^||ae(7o(G^)||oo <oo 



Now we recall ( 2.21 ), estimate the absolute value of the Poisson brackets by 
the products of the supremum norms of the derivatives, employ Lemma |2.8| and 
conclude that the contribution of the Poisson brackets involves {G{B))™, where 
C{B) depends on || V^cro(i?)||co- Thus for a certain G3(i?) < 00 and any L eN 

Tslz^Ksf*^)) -Tgiz^KB) <C3,L(S)-m2.(G3(B))".M-^, M > 1, 
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where Cz[B) depends on || V^cro(-B)||oo and || sub(i?)||oo- The factor m? is du e 
to the fact that there are m terms m the expansion of the Poisson bracket in ( 2.43 ). 

3. Contribution of ao[B^^ ■ ■ ■ B^^) 
Let B e 'I'"(Af) (we redenote G by B). Our goal here is to calculate 

(3.1) (^™(^))" / MB^)dxd^ 

xiH l-x„=0 -JS'M 

for n — 1,2. Assume for the moment that n is any natural number. Recall 
that Sm denotes the set of all permutations r of the numbers 1, • ■ ■ ,m. Denote 
Jc-r := (><ri, • • ■ , ^Trr^)- Note that both the second factor in (3.1) and the domain of 
summation are symmetric in Therefore 



-. m— 1 „ 



>ti^ h>tm=0 reS„ p=l 

Now an application of the gHD (Theorem 8.2) gives 

m-l min(p,n) 



{~{k,{Hr))-y^ (-(fc,(3?0)-)'^ 



TO! ' ^ ' ^ ' V. 

KiA \-Km=0 r&Sm P=i k^, --.kj>l ll,- -.!j>l 

k-^^ ^k;=p \-l^=n 



n 



ki 



o'oiB^) dxd^, 



S'M 



where the notation (S.l) has been used. We make the sum ^ X^tgS ^^'^ outmost 
one and then drop it, because it is the number of >f's that matters, not their indices 
(symmetricity of the last factor and of the domain of summation again). We get 



m — 1 min(p,n) 



(3.2) 



-:Kj-n—0 p—1 J — 1 



,fcj >i ii .<j >i 

ki+--- + kj=p ij-|-...-|-i^=„ 



n 

h, - ■ ■ ,L 



(-(%(3?))_)'^- 



o'o{Bj^) dxd^. 



kl i^j J S'M 

We have m — 1 independent summation variables Jc. Let us make the change of 
variables 

:^ ki{J() ^ >ii -\ h 

vi ■■= k2{'>c) = Xkt+i H h Xki+k2 



Uj kj{K) = >Cki+...+kj-i+i H 1- >Cki+...+kj_i+kj 

for j — 1, ■ ■ • , min(p, n). We also interchange the summations over p and j, intro- 
duce the new variable fcj+i := m — p, and perform the to — j — 1 unconditional 
summations in ( |3.2| ) . We need the following important fact [ 
which states that the convolution on the Fourier transform side corresponds to rais- 
ing S to a power on the original side. 



G02, after Lemma 1.2] 
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Lemma 3.1. For any B e 5'°(M) and arbitrary j G N and v ^% one has 



>^iH V>tj=v 



Notice also that by Egorov's theorem ( pj] ) and because (7n{B^) = {an{B)Y , 



Then (3.2) becomes 

min(n,m— 1) 



1 ' t,^ ... 7, ■ 7 , ... 1 - - 



3 = i 



(3.3) 



Z^l , ■ ■ ■ ,Z^j i 1 , ■ ■ ■ 7^ J > 1 



27r 



n 

h, - ■ ■ Jj 



2tt 



X / dxd£_ 



2tt p2tt 



t{iyiti + ---+^jtj-{iyi + --- + iyj)tj + i) ^ ^ ^ , d,tj dtj+i 

27r 27r 27r 



"'0 

xF,[z™](fe[,^...,6*^5*^+M 

where we have mtroduced a hnear [j + l)-map i^j+i, whose action on f{z) = z™, 
m > J + 1, is prescribed by 



(3.4) Fj+i[z"']{xi,--- ,Xj,Xj+i) 



and Fj+i[z™] := for m = 1, • • • , j. 



fel >■■■ .kj ,fcj+i >i 

fclH l-fcj+fcj^i=Ji 



^ l2_ + : 

fci ■■■ ^- ""j+i 



Remark 3.1. It is important that for any fixed n the sum over j in ( |3.3D terminates 
at n, no matter how large m is. 

We need the following fact pO% after (1.3)]. 



Lemma 3.2. Let M be a Z oil manifold and A definedby (1.1). Assume (7i{A){x, — 
cri{A){x, for all {x,^) G T*M. Let f,g e C'^{S*M), v and recall the no- 
tation ( |l.2| ). Then 



S'M 



f„g-„dxd^^ / f_i,g^dxd^. 



S'M 



We need to compute (3.3) only for n — 1, 2. The computation for n = 1 has been 
carried out in [ G02 |. We repeat it here for the sake of completeness. For n = 1 
(B.2) becomes 



2tt p2tt 



E(-(^i)-) / dxd^ r I 

JS'M Jo Jo 



/ . L y°o ) ~~ 



(3.5) 



ki ,k2>l 
ki+k2 = m. 



ki 



27r 2tt 



S'M 



fel >fe2 >1 
fel+fe2=m 



This can be rewritten as 
(3.6) 



E(-M-) / dxd^ E - c^o^')-^ 

Js'M ^.^1 \3 



(&0^)-.i 
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which by Lemma 3.2 equals 



(3.7) 



m-i ^ ^ , J 



S*M 



where the last equality is obtained after the change of summation index j m — j . 
Now we sum (^) and (3/7) and use i + -^^^ = and Lemma 3^ to conclude 

that (3.5) equals 



m 
' 2" 



„ m — ] 

^(-(^)_) / dxdi 

JS'M .1 



fe=i 



i=i 

27r /.27r 



dxd^ 



"'0 



2^ 2^' 



where following LRS[ | we have noticed that W2 defined by (1.5) satisfies 

m — l j m~j 



J TO-J 



Now for the case when n = 2 in (3.3). In this case j can be either 1 or 2. In the 
first case we repeat the above argument, the corresponding contribution becomes 

'■^'^ '■^'^ ,ti ,i2N dti dt2 



E^^ 

fe=i 



S'M Jo Jo 



*fc(ti-t2)p^2[z™](&^\5*^) 



27r 27r 



(3.8) 



For j — 2, ( |3.3| ) becomes 

Y.{-{^l)-){-{^2)-) I dxdS, 

where by (3.4) 

(3.9) f^3[2n(2^i,2;2,X3) := 



m > 3, 



and F3[z™] := for m = 1,2. Let us redenote W3 := F3. The formula (|l.8D for 
W^si/]? / G is justified in Secti on |7[ Now the last detail, the following can be 
proved in the same way as Lemma \i.2\ For any f,g,h& C°°{S*M) and A, /i, G Z 



(3.10) 



/a '^^'^C = / /-A ^'-1^ c'sjd^. 



In view of (3.10) we rewrite ( |3.8| ) as 

r.27r /.27r p2tt 



kl 



dxd^ 



IS'M 



Jo Jo 



/TT ZTT /TT 
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4. Contribution of the symmetric part of suh[B^-^ ■ ■ •i?>c,„) 

Let us choose any B e 5'°(Af). The goal of this section is to calculate for an 
arbitrary m e N, to > 2, the sum 



(4.1) 



V M„(>?) / ( 



dxd^ 



>€l-\ \->C,„=0 



^ <Jo{B^, ■ ■ ■ B^^_J sub (B^^) ao{B^^^^ ■ ■ ■ B^^), 



r=l 



As in Section ^ the key observation is that both the second factor in (4.1) and 
the domain of summation are symmetric in >?. We permute all the >r's in the first 
factor in (3.1), make use of Theorem 8.2 for the power n — 1, that is the classic 
HD, 

^{Xti H 1- >^T, )- 



J 



and drop the summation over Sm, because, as in Section |^, it is the number of x's 
in a particular sum that matters, not their indices (again, the domain of summation 
and the second factor are still symmetric). After that we get 



d-d^ EE E 

j—l r— 1 >fi-) \-:i^jji—0 



-(xi H h Xj)- 



X cro(-S^i • • • S>c,._i) sub (B^,,) cro(S^^+i • • • S^,„). 

We split the above sum into two: for r = 1, • • • , j and for r = j + 1, ■ ■ ■ ,m. In the 
first case we set 



Then 



><1 + • • • + >ir-\ + Xr+l + • • • + = ^, + l + 



We carry out the independent summations and recall Egorov's theorem (2.1) and 
Lemma 3.1 to write the sum over ?' = 1, • • • , j as 



(4 



,2) / dxde E -EE(-(^)-Vo((^'"')-M)sub(B^)ao((i?'"-^)-.) 

Now let us write -.^ (7o{H){e\x, ^)) and /i'^^ := sub(iJ)(e*(a;, ^)), He ^"(Af). 
We have 



2-R 



'ib'oy~'K^^-^-suh{iBr),) 



2tt j 



Then (4.2) becomes 



(4.3) 



/ dxd^Y -E(-M-)«'^b((i?^),)ao((i?"--'")-.) 

JS'M J ^ 
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Analogously, in the case r = j + 1, ■ ■ ■ , m we set 

v := >Ci + ■ ■ ■ + >Cj , —fi ~ V :— x^.. 

Then 

+ 1 H \- Xr-l + Xr + 1 H ^ ^l, >C.j + i H h >fm = -V. 

Then the sum for r = j + 1, • • • , m becomes 

m — 1 m 

„ m—1 

which in view of Lemma ^.^l equals 

(4.4) / dxd^ E - sub ((S"-^).) ^o((i?^)-.) 
We make now in (4.4) a change of index j m — j and get 

^ m—1 

(4.5) / dxdC E(-(^i)-) E ^ sub((B^).)fTo((i?'"-^")-.) 



We sum (gj) with gj), use i + 



j{m-j) ' refer to Lemma 3^, and replace 
j m — j again to conclude that (4.1) equals 



(4.6) 



S*M 



rn — 1 

E^+E 



sub 



It remains to notice that ( [4.q ) in view of (2^) equals 

/ dxdi E W / ^2 [;^"] (6^^ , hl^ ) 



dxd^. 



dti dt2 
277 2?:' 



5. Contribution of the non-symmetric part of sub(i3^j • • • -Bj^„J 

Let us choose any B G \1/''(M). The goal of this section is to calculate for an 
arbitrary m € N, m > 2, the expression 
(5.1) 

P m 

>^l + -" + >t™=0 JS'M l<k<l<m P=i 

The domain of summation is symmetric in 7c. However the integrated over S*M 
sum is generally speaking not symmetric. Each of the terms in the sum in 
possesses a partial symmetrj^^ 

Depending on k and / in (O) there are three possible cases. In the first case all 
the indices except for those in the Poisson bracket form one continuous block. In 
this case it is convenient to rewrite Mm(3<) as follows. Note that for a, 6 G M 



(5.2) 



min(a, b) = a — {b — a)_, 
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let p e N, p > 2, /ii, • • • , /ip g Z and set Sj :— + ■ ■ ■ + fij, j = 1, • • • ,p. Then for 

any j = 1, • • • ,p- 1 

min(0, si, • • • , sj, Sj+i, ■ ■ ■ ,Sp) 

= min( mm(0,si,--- , Sj), min(sj , s^+i, • • • ,s,„)) 

= min ( min(0, si, • • • , sj), Sj + min(0, Sj+i ~ Sj, - ■ ■ , Sp — Sj) ). 

Therefore it follows from (5.2) 

Aip(^i, • • • , /ip) = Mj{^i, ■■■ , ^ij) 
(5.3) , . 

- (Ml + ^Mj - Mj{fii, ■ ■ ■ ,/ij) +Mp_j(^j+i, • • • ,/ip)j_- 

The three subcases of the first case are 

A; = 1, ^ = 2; k ~ 1, I — m; k = m — 1, I ~ m. 

In these three subcases using ( p.3| ) we rewrite Mm{xi, • • • , >?,„), respectively, as 

M2(>!ri,>f:2) - (>f:i + >f2 - M2ixi,>{2) + Mm-2(>f3, ■ ■ • ,Xm))_, 

and 

-(>!ri)_ - (xi)+ + Mm_2(>*2, ■ • • , Xm-l) 

- { >C2-\ h ><,„_1 - M„_2(>f2, • • • , >^m-l) + >fm ) _ ) , 

and 

Mm_2(>!ri, ■ • • ,>irm-2) - ( ^^1 H h ^<m-2 

- M„_2(>*1, ■ • • , Xrn-2) + M2i>Cm~l,Xrn))_- 

In the second case all the indices but the two from the Poisson bracket form two 
continuous blocks. In this case we rewrite Mm(>?) as follows. Note that for p e N, 
p > 3, ■ ■ ■ , fip and si , • • • , Sp as above, and any j = 1, ■ ■ ■ ,p — 2 

min(0, si,-- - , Sj, Sj+i, • • • , Sp) 

= min ( min(0, Si, • • • , Sj), min(sj+i, • • • , Sp) ) 

= min ( min(0, si, • • • , Sj), sj+i + min(0, s^+i — sj, - ■ ■ ,Sp — Sj) ) 



which together with (5.2) implies 

(5 4) ^^P*-^^' " ' ' " Mj{>ci, ■■■ , fij) 

- (mi H \- f^] - Mjifii, - ■ ■ fij+i + Mp_j_i(^j+2, • • • ,//p))_. 

There are three subc ases: first fc = l,^ = r + l,2<r<r7i — 2, in which A/,„(>?) 
after first using (^.3|) and then (|]J) becomes 



-(>^l)- - ( + Mr-l(>f2, • • • ,><r) - (><2 H h ^f,, 

- Afr-l(>?2, ■ ■ ■ ,Xr) + Xr+1 + M,„_r-1 (>fr+2 , ' ' ' , >fm) )_ ) , 

second: A: = r + l,Z = r + 2, l<r<m — 3, in which A/m(>5) after first using ( ^.4| 
and then (|]^) becomes 

Mri^i, • • • , ><:,.) — (^Xi + ■ ■ ■ + Xr — Mr{>Cl, • • • , >fr) + Xr+l 

- {Xr+2)- - ( {}<r+2)+ + Mm-r-2iXr+3, ' ' ' , >fm) )_ ) , 
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third: k = r + 1,1 = m, l<r<m — 3, in which Mm{x) after first using first (5.4) 
and then (5.3) becomes 



Mr{Ki, • • • , Hr) — (^Xl + ■ ■ ■ + Xr — Mr{>Cl, • • • , >Cr) + Xr+l 
+ Mm-r-2{Xr+2, ' ' ' , >*m-l) — ( >*r+2 + ' ' ' + >*m-l 
— Mm_r-2(><r+2, • • • ; >^m-l) + ^^m )_ ) ■ 

Finally in the third case the indices from the Poisson bracket are taken to be 
>f,._i_i and Xs+i, where 

l<r, r + 2<s, s + 2<77i, 
and there are three continuous blocks 1, • • • ,r and r + 2, • • • , s and s + 2, • • • , m. 



In that case Mm{x) after using (5.4) twice becomes 

Mr{xi , • • • , — + • ■ ■ + ><r — Mr{xi , • • • , >fr) + ^i'r 



+ 1 



(5.5) + Ms^r-l{Xr+2, ' ' ' ,><s) ^ ( ><r+2 H h >fs 

— Ms_r-l(>*r+2, • • • , Xs) + Xs+1 + M„i^s-l{Xs+2, ' ' ' , ^m) , 



We make the computation for the third case, the first and the second are treated 
in the same way. We use the a convenient reformulation of the original form of 



the BSt (Theorem B.l), and one property of the j-maps j G N, see Lemma 5.1 
below. We rewrite the corresponding to the third case part of (|5.lD as 



m — 3 m — 2 



(5.6) 



p III/ U I I L JLi 

■ ■^^ =n J S' M s—r+2 



>tiH h>^„.=0 



X o'o(S>fi • • • B^^,) <To{B^^^2 ■ ■ ■ B^J o-o(i?>c^^2 ■ ■ ■ B^, 



where k — r + l, I ~ s + 1, and each of the three products under the (Jq sign 
contains at least one factor (therefore this expression is non-zero only for m > 5). 
Now we observe that both domain of summation and each of the three products are 
symmetric if we interchange the indices xi, ■ ■ ■ x^, >i'r+2i ■ ■ ■ j and >^s+2J ■ ■ ■ i 
separately, that is within each of these three sets. Then their sums do not change, 
and we can therefore consider the representation as a function of 

Mr{xi,- ■ ■ ,Xr), Ms-r-l{Xr+2,- ■ ■ ,Xs), Mm~s-l{Xs+2, ' ' ' , Xm) 



.1 



and 



only. We interchange the indices within the three groups, use Theorem 
take in account Remark 8.2. After that because of the symmetricity of domain of 
summation, and of the principal symbols of the three products we may, and will, 
drop the coefficients l/{r\), l/((s — r — 1)!) and l/((m — s — 1)!) and the summations 
over the permutations over the three groups of indices (just as in Section 0). We 
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conclude that (5^) becomes 
(5.7) 



m — 3 m — 2 r ^ ^ s — r—1^ 



yc-^^ r— 1 s— r4-2a— 1 



3\^ \-3a=-r 



(3=1 



l-lH |-fc^=S-T- 



ki---kp 



m — s — 1 



7= 



ilH K-, = m-s-l 

(jl (>?))_ (ia (>?))--( (jl (><))+ H H (ja + 

- (>?))_ (fc/3(>5))- - ((/ci (>?))+ + •• • + [kp[Tc))+ + m:s+i 



X o'o(5>fi ■ • • B^^) (To(i3j^^^2 • • • i?>cj (^oiB^^^^ ■ ■ ■ B>t,„)- 

Here for all possible values of the indices eac/i of the summation variables >fi , ■ • ■ , 
is involved in some of the jiijc), ■ ■ ■ , but only once. We rewrite now the sum over 
r and s as a sum over three summation variables a,b,c> 1, which are the lengths 
of the three continuous blocks, with the condition a + b + c — m — 2. Note also 
that it does not matter anymore on which place >Cr+i and >fs+i stands, it is only 
important that the same letter is used in the non-positive valued function coming 
from the square bracket in ( 5/7 ) . Let us make two changes of variables 



Xl :— Xr+l, X2 ■= Xr+1 

and for all a,h,c> 1 with a + 6 + c = 77i — 2 

/ii (>?),••• :=jQ(3<), l<a<a 
i^i := fci(3?), • • • , := fc/3(>?), 1 < /3 < 6 
pi := • • • := 1 < 7 < c. 



Then the square bracket in (p.7|) becomes exactly the defined by (1.12) non-positive 



valued function 



^(3) 



Now we carry out the ((a — a) + (6 — /3) -t- (c — 7)) free 



summations over the rest of >ir's, just as in Section ||. We get 



(5.8) 



a b c ^ - ^ 



E 



^l + '*2+MlH + 

+ ■'! + ■ ■■+i'3+PlH hP-,=0 



{CTo(-B>ri),cro(-B>c2)} C^O ( E 

\ ii, - - ,JQ>1 

JlH hja=a 



ii • • • 



LOWER ORDER TERMS IN SZEGO TYPE LIMIT THEOREMS ON ZOLL MANIFOLDS 33 



X do 



E 



E 



Pi 



fclH l-fc/3 = '5 



'l-l h'7=<: 



• • • L 



Now by the definition of the Fourier coefficient, by Egorov's theorem and the fact 
that ao{B^) = {a(){B)y , j € N, we obtain, just as in Section ||, 



dxd£_ 



S'M 



LL U X 1 1 

E EEE^flTry 

,,b,c>l a=ll3=li=l ' ^' '■ 



E 



(5.9) 



>£l + >t2+MlH + 

{ao(S^J,ao(S^J} 

gi(;iiriH h^iarc+i^isiH hi^f3Sa+pitiH \-p-,ty) 



dt-y 

2^ 



X [z-^] , . . . , 65° ) $^ [z"] . . . , 6^'^^ ) ci>^ [z^] (, , . 



6*^ 



We use the fact that <i>Q[2;''] — for a > a to extend the summation over each of 
the variables a, (3, 7 to the whole of N, and make these summations the outer ones. 
Finally we use the following important property of the j-maps $j, j € N, whose 
proof immediately follows from the definition of $j and is omitted. 

Lemma 5.1. Let choose an arbitrary p G N. For any fixed real or complex Xj, yk, 
zi, j, k,l = I, ■ ■ ■ ,p — 2 the following identity holds 



VVV— — — 

Q,i «i ^1 

1 Q=l/3=17=1 ' ^' 



a,b,c> 
a+b+c=p 



(5.10) 



X $„[Z°](X1,--- ,Xa)^f3[z'']iyir-- ,2//3)$7[2;T(zi,--- ,Zy) 



00 00 00 



E E E Q,i /31 ^1 

a=l/3=l7=l ■ '■ 



X $a+/3+^[zP](xi, • • • ,Xa,yi,--- ,yp,zi, ■ ■ ■ ^z-r)- 

Note that the sums terminate for a + /3 + 7 > p, and a, /3, 7 > 1. Therefore only 
p o f the variables xi, • • • , Zp-2 are present at the each term on the right-hand side 
in ( |5l0|) . 

The formula ( |l.l5| ) now follows from ( |5.9| ) and Lemma |5.l| with p = m — 2. It is 
important that we have reduced the number of j-maps from three to one, and that 
the j-map ^j[f] is linear in / for all j G N. Therefore (1.15) holds for any f G Ai- 
We divide by z^, because for a monomial /(z) — z™, m > 5, the j-map should be 
evaluated at z™~^. We subtract the fourth degree Taylor polynomial T4[f] of / 
about t = because the term under consideration is absent for all polynomials / 
of degree < 4. 

The formulas (1.13) and (1.14) are proved analogously. Because the correspond- 
ing expressions appear for the monomials of degree at least 3 and 4, respectively, 
we subtract T2[f] and T^lf], respectively. 
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6. Proof of Proposition 1.4 



We will need the following statement (see for instance | GrRy ). 
Lemma 6.1. Let n,TO e N. Then the following holds, as n oo, 
" n"^^^ 1 _ TO .„ 



^ m + 1 ' 2^ ■ 12" 

A;— 1 

n ^ 

^ fc"^ = logn + 7 + -n"^ + 0(n"^) 



0{n" 



k=l 



n 



C{2)- -+0{n-^) 



= C(to) + 0(71-™+^), TO>3. 



fc=i 



Note that 
(6.1) 



Tr(P„G) =^Tr(7rfcG), 7i e N. 
fc=i 



Assume first that d > A. Then by (|2J) 

Tr(7rfcG)-fc''-ii?o(G)-fc''-^i?i(G)-fc'*-^i?2(G) < ci(G) fc'^"'^, /c G N, 



and so in view of (3.1), for any n e N, 



Tr(F„G)-n'^- -i?o(G)-n'*-i 



ii?o(G) + ^i?i(G) 



+ C3(i?o(G),i?i(G),i?2(G),d) < c2(G,d)n'^-3^ 



(6.2) 



For d = 3 there also appears a term with logn in the left-hand side of ( |6.2[ ). This 

proves part (iii). 

For d = 1, 2 when we sum over k — 1, • ■ • ,n in (|2.9|) there is a subtle point, 
nam ely, the constant coefficient in (O), as n ^ cx). The terms of all orders in 
( 2^), and also the possible rapidly decaying term, will contribute to it. Assume 
(1.17). Then the following series is absolutely concergent 



(6.3) 



efc(G) Tr(7rfeG) - ^ k^-'~' Ri{G), keN. 



1=0 



Furthermore, for any N eN and all fc e N by ( |2.9| ) 

N+d~2 oo 

|efc(G)| < |Tr(7rfcG)- ^ k^-'-' Ri{G)\ + \ ^ fc"^-!-' i?,(G)| 



(6.4) 



l=N+d-l 



< cn{G) fc"^ + fc-^ |i?z(G)| < cXG) • fc- 



i=0 
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in view of ( p7| ). Note that C{G) defined by ( p^ ) equals 

OO 

(6.5) C(G)=^efe(G), 



fe=i 



the series being absolutely convergent. Now summing over A; = 1, 
obtain 



, 72 in (|6.3| ) we 



+00 



(6.6) 



Tr(P„G) = J2 ^'(G) E fc'"'"' + E ^^(G)- 



i=0 



fc=i 



fc=i 



Because efc(G) decays rapidly ( |6.4|), a s fc — > 00, and by ( |6.5| ), X]fc=i £fc(G') converges 
rapidly to G(G). From this and (6^) we can obtain an asymptotic s of Tr(P„G) up 
to any negative order. In particular, taking into account Lemma |6.l| we prove (i) 
and (ii). 

7. A FORMULA FOR $j [f], f £ Ai, j eN 

We find a formula for $j [/] , / G Ai , in terms of an auxiliary linear J-map $ j , 



j € N, see ([7.2[). The latter acts on monomials as follows 



(7.1) 



^.■[^'"Ka,--- ,C,):= E ^---^^ m = 0,l,2,..-. 



il, -- ,ij>0 



which is the complete symmetric function of degree m evaluated at the point 
(Cij • ■ ■ J Cjj 0, • • • ). To write a formula for ^j[f], f G Ai, we use an idea sug- 
gested by Kurt Johansson. It uses a Cauchy integral representation of (7J_) via the 
generating function. Namely, by the identity 0, (1.2.5)] for {£,!,■■■ i Cji Oi ' ■ ' ) 



E*.'-[^"K6,--- ,0 



n 

k=l 



C-6 



After a multiplication by C™ ^ a-nd an integration over a contour 7 C 
circumferences the points 0,^i, • • • we single out <tj[z™'] and obtain 

1 f 



which 



27ri 



(c-a)---(c-c,) 



Therefore for / G Aq (an analytic on C function which might have a constant term) 
we can define 



1 

27ri 



(c-a)---(c-6) 

The j-map $j is now defined for any / G Ai by 



(7.2) <i>,[/(z)](xi,. 



*.-[^"V(^)](6 



fiOdC 



Now we take ( |7.2| ) as a definition of $j, then (1.9) with m > j holds, and it only 
remains to prove that $j vanishes on the set of polynomials of degree j — 1 with 
no constant term. Note that the integrand in (7.2) equals 



(7.3) 



1 

27ri 



rV(c) 



(C - Ci) • • • (c - ij) 
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The integral does not depend on the contour of integration, if only all ^i, • • • , S,j 
are inside it. Let 7 = |C : |C| = R}, R — > co. If f{z) 



1, 



• , j — 1, then 
-y 0, as 



the absolute value o f (|7.3| ) is estimated by {2tt)-^R^-'^ ■ 2ttR/R^ = R 

i? — » 00. Therefore (7^) must be 0. 

If / G .4i then z^^f{z) is analytic. Therefore (7.2) equals the sum of j residues 
at the points ^1, • • • For instance, for j — 2 and any / e Ai 



*2[/](xi,a;2) = 



4'V(6)-g2"V(6) 



Remark 7.1. Another way to construct j € N, is by induction on j. In that 
case one uses a simple formula 



(7.4) 



,r+l 



„r+l 



= u^-y*, u,v G 



p,q>0 
p+q=r- 



for r = 0, 1, 2, • • • . This derivation is however longer than the above argument. The 
formula (7.4) was used in |LRS| for the computation of the 2-map W2 defined in 
(1.5). In that case the induction is not needed. 

Finally we find a formula for M^3[/], f E Ai, where the action W3[z™], to e N, 
is given by (3.9). In view of (7.1) using the integration as in (7.2) and moving out 
X3 we obtain 



W3[f{z)]{xi,X2,X3) = X3 



$3[^~V(^)](a, 6, 2:3)^6^6- 



Now evaluating the three residues in (7^) for j = 3 we obtain (L 



8. Formulation of the gHD, and of a version of the BSt 

In this section we state the gHD and a convenient for our purposes version of 
the BSt, see ]Gi3{ for details. Let to e N and xi, . . . , >Cm G M. Recall the notation 
( p.4| ), and for each permutation r e Sm write 

Fix any j — 1, • ■ • , m. Recall [Coh| that a partition of to is a way to write to as 
a sum of natural numbers with no respect to the order of the summands. These 
summands are called parts, and their values are called lengths. Note that each of 
the parts has length at least one. For an arbitrary partition of to into j parts 

ki > 1, ■ ■ ■ , kj > 1, ki + ■ ■ ■ + kj = m, 

we introduce the notation 



(8.1) 



ki{Kr) := Xt 
k2{jCr) 



Each of ki{7c^), I — 1, ■ ■ ■ , j , is a sum of fc; permuted variables out of Xr^, • ■ • , tir^ 
so that each of the permuted variables enters exactly one sum. Note also that 
because ki + ■ ■ ■ + kj = m one has 

kl (TCt) + ■ ■ ■ + k-j iJir) = Xti + ■ ■ ■ + >*r„ 
= Xi H h Xm- 
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For a S R, denote — (a)_ := min(0, a) and (a)+ := max(0,a). We state now 
a combinatorial formula, called CF, which is an equivalent form of the BSt |^ 
Theorem 2.2], and does not involve any advanced combinatorial coefficients. This 
formula is very suitable for the calculation of sums which arise in Secti on [s] and in 
computations of the joint distributions for random walks as in Theorem 1.7. Recall 



the notation (2.4) 



Theorem 8.1 (CF: an equivalent version of the BSt). For any m G N, arbitrary 
, Xm G M, and any real- or complex-valued function f defined on the left 
half-axis, the following holds 



.2) " 1 



re5,„ j=l 



fci , -■ ,fcj >l 

fclH hfe,-=m 



f{-{h{Kr))^ (fc,-(^r))-) 



This holds because the sets of the values of the arguments of / on the left- and on 
the right-hand side in (B.2) contain the same numbers with the same multiplicities, 
by the BSt, the rest being just an account of the number of conjugacy classes in 
Srm see Section 3 and especially the proof of Lemma 3.3 in |Gi3|. 

For a monomial f{z) — z", n G N, a further calculation in CF can be carried out, 
of one subtracts the value of / at the "previous" maximum. The right-hand side 
of the resulting formula has a multiplicative, and not additive, as in the CF, form, 
which is important for the calculation of the sums of convolution type in Section ^ 

Theorem 8.2 (Generalized Hunt-Dyson formula). For any power n G N, an arbi- 
trary number of variables m G N and for arbitrary fixed values of the real variables 
xi, - ■ ■ , Km one has 

J2 [(M™(>?r))"-(M™-l(^r))"] 
min(m,n) _^ 

~ E E E 

res™ J=l 



(8.3) 



fcl .fsj >1 !l , -■ j >1 



E 

•■•+!_ 

kj 



n 

li, ■ ■ ■ ,lj 



{-{k,{->ir))-r 

kl ■ 

Remark 8.1. In the case n ~ 1 the identity ( |8.3| ) becomes the usual Hunt-Dyson 
formula @, after (4.8)] 



(8.4) 



J2 [ Mm (>?r ) - Mm-lC^r)] 



res,. 



{m-l)\{^{Kl + --- + Km)-). 



Remark 8.2. It is important that Theorem B.l and B.2 holds if the symmetric group 
Sm is replaced with a larger symmetric group of indices. One does however need 
the summation over the whole group Sm ■ The usual formula HD holds for the sums 
over the cyclic subgroup Cm, as well. 
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Remark 8.3. The usual formula HD of course can be viewed as a corollary of the 



Bohnenblust-Spitzer identity. However one needs the formula (B.2) for all n £ N 



to reprove the BSt as is done in |Gi3, Section 4] 



Remark 8.4. The BSt, the gHD, and the HD hold if the minima are replaced with 
the maxima, and the negative parts — (•)- replaced with (•)+. 
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